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Mirror dark matter is a dissipative and self-interacting multiparticle dark matter 
candidate which can explain the DAMA, CoGeNT and CRESST-II direct detec- 
tion experiments. This explanation requires photon-mirror photon kinetic mixing 
of strength e ~ 10~ 9 . Mirror dark matter with such kinetic mixing can potentially 
leave distinctive signatures on the CMB anisotropy spectrum. We show that the 
most important effect of kinetic mixing on the CMB anisotropies is the suppression 
of the height of the third and higher odd peaks. If e ~ 10~ 9 then this feature can be 
observed by the PLANCK mission in the near future. 
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A large variety of observations have lead to a simple picture of the Universe. In 
a nutshell, we live in a spatially flat, expanding Universe consisting of dark energy 
(~ 70%), non-baryonic dark matter (~ 26%) and ordinary baryons (~ 4%). In the 
last two decades, detailed observations of the cosmic microwave background (CMB) 
by COBE[l], WMAP[2], SPT[3] and many other missions have provided important 
tests of this basic picture. A key question concerns the identity of non-baryonic 
dark matter. Although it is popular to assume that dark matter consists of a single 
species of weakly interacting massive particles (standard cold dark matter model), 
in actuality the particle physics underlying the non-baryonic dark matter in the 
Universe is currently unknown, as is the physics responsible for the dark energy. 

One thing we do know, though, is that the standard model has been very suc- 
cessful in describing the interactions of the ordinary particles. In fact, it is possible 
that such a structure might also be responsible for the non-baryonic dark matter in 
the Universe as well. That is, dark matter might consist of a hidden (mirror) sector 
with particles and interactions exactly isomorphic to the ordinary ones[4, 5] (for a 
review, see e.g. [6]). Provided that initially the mirror sector temperature is much 
less than in the ordinary sector, i.e. T' <C T in the early Universe, such a scenario 
can explain the large scale structure of the Universe in a way completely analogous 
to standard cold dark matter [7, 8]. 

On much smaller scales, mirror dark matter is radically different to standard cold 
dark matter. It is self-interacting, dissipative and mult i- component 2 . These proper- 
ties might help explain some puzzling aspects of dark matter on small scales, such 
as the inferred cored central density profiles in galaxies [c.f. ref.[ll]]. At the current 
epoch dark matter needs to be roughly spherically distributed in spiral galaxies to be 
consistent with various observations 3 . However, within galaxies mirror dark matter 
would be expected to collapse into a disk, analogous to the way in which ordinary 
matter collapses, unless a significant heat source exists. It turns out that ordinary 
supernovae can potentially supply the required energy provided that photon-mirror 
photon kinetic mixing[12, 13], 

r . — - f^F' (i) 

of strength e ~ 10~ 9 exists[14]. [In the above equation, F^ u (F' ) is the ordinary (mir- 
ror) photon field strength tensor]. That is, the ordinary and dark matter components 

2 There are potential limits on self interactions of dark matter from observations of the Bullet 
cluster [9]. These observations can set stringent limits on dark matter self interactions provided 
that the bulk of the dark matter particles are distributed throughout the cluster and not bound to 
individual galaxies. However mirror dark matter is dissipative and in clusters (or at least in some 
of them) the bulk of the dark matter particles might be confined in galactic halos [c.f. ref.[10]]. 
Under this assumption mirror dark matter is consistent with Bullet cluster observations. 

3 Deviations from perfectly spherical halos are also required and might be due to various sources 
including, e.g. partial collapse of the halo due to dissipative processes, mirror magnetic fields, 
asymmetric heating from ordinary supernovae which are distributed in the disk (i.e. they are not 
spherically distributed), possible existence of a dark disk subcomponent etc. 
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of galaxies might be in a kind of dynamical equilibrium where the energy supplied to 
the halo by ordinary supernovae balances the energy lost to the halo due to radiative 
cooling. Since the former is related to the galactic luminosity and the latter the dark 
matter density, that is, v rot , it has been speculated [15] that this might potentially 
explain puzzling regularities on small scales, such as the Tully-Fisher relation[16]. 

Photon-mirror photon kinetic mixing of strength e ~ 10~ 9 is also implicated[17] 
by the positive results of the direct detection experiments, DAMA[18], CoGeNT[19] 
and CRESST-II[20]. It has been shown in ref. [21] that these experiments can be 
explained by the interactions of a halo Fe' component with ey^v ~ 2 x 10~ 10 . Here 
£, Fe > is the abundance by mass of the halo Fe' component (at the Earth's location) 
normalized to 0.3 GeV/cm 3 . Naturally, it is very difficult to predict £p e / with any 
certainty, but of course we expect £p e / ~ 1, and thus e ~ 2 x KT 10 . 

Kinetic mixing can also have important implications for cosmology. Successful 
big bang nucleosynthesis (BBN) and large scale structure (LSS) require the initial 
condition py <C p 7 and ny ~ 5n&. How such initial conditions might have arisen has 
been discussed in the literature [22]. However, if photon-mirror photon kinetic mixing 
exists then this will generate entropy in the mirror sector via the process ee — » e'e' 
when T 7 ~ m e [23]. That is T 7 will be generated even if we start with T 7 <C T 7 . In 
fact, it has been shown that the asymptotic value of the ratio: T 7 /T 7 , which we here 
define as x, is given by [24, 25] 



The value x ~ 0.3 is close to the limit estimated from the matter power spectrum i.e. 
successful large scale structure [8]. We will estimate that the upper bound on x from 
such considerations is conservatively around x ~ 0.5. Non-zero T 7 /T 7 will also lead 
to important effects for the CMB as previously discussed in ref. [8] (see also ref. [7]). 
In view of the forthcoming results from the PLANCK mission it is pertinent to 
examine thoroughly the possible effects that kinetic mixing will induce for the CMB. 
This is the purpose of this paper. In fact, we will show that mirror dark matter can 
potentially leave a distinctive imprint on the tail of the CMB anisotropy spectrum. 
Our most important observation is that the height of the third and higher odd peaks 
can be suppressed. This should be observable by PLANCK provided e ~ 10~ 9 . 

To summarize, we assume a mirror sector exactly isomorphic to the ordinary one, 
except with initial conditions <C T 7 . Ordinary and mirror particles can interact 
via photon - mirror photon kinetic mixing, which can excite the mirror degrees of 
freedom in the early Universe. In particular the process: ee — > e'e' generates the 
mirror particles until the ee have annihilated at T 7 ~ m e , the final T 7 /T 7 value given 
in Eq.(2). In fact, most of the mirror entropy generation occurs after the neutrinos 
have decoupled. One effect of this is to induce a slight cooling of the ordinary photons 
relative to the ordinary neutrinos. The net effect is that there is additional neutrino 
energy density and also an additional relativistic component comprised of mirror 




(2) 
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photons. These two additional components to the relativistic energy density can be 
parameterized in terms of extra neutrino degrees of freedom: [25] 



Here, the temperatures are evaluated at the time when photon decoupling occurs, 
i.e. when T 7 = T dec ps 0.26 eV. Using the result from Eq.(2), together with the usual 
Tj,/T 7 relation, we have 



Also, numerical work[25] has found that 5N^ ff (e)[CMB] « 0.8N% ff (e)[CMB]. This 
additional energy density can directly affect the predicted CMB anisotropics. In 
fact, it is known that additional relativistic energy density can dampen the CMB 
tail[26, 27]. However, there is another important effect for the CMB. If dark matter 
consists of mirror particles, then they experience significant pressure prior to mirror 
photon decoupling. If T 7 < T 7 then this epoch occurs prior to the familiar hydrogen 
recombination. One can anticipate that the small scale inhomogeneities in the mirror 
matter density should be suppressed, since the Fourier modes which enter the horizon 
before the time of mirror hydrogen recombination undergo acoustic oscillations due 
to the pressure of the mirror baryon-photon fluid 4 . In other words, we expect a 
suppression of power on small scales when compared with standard non-interacting 
cold dark matter. The previous study[8] has indeed observed this effect on the matter 
power spectrum. This suppression of power on small scales will also influence the 
CMB spectrum, and one would anticipate that this might also dampen the CMB 
anisotropics at high multipoles. This effect, is of course, in addition to the effect 
of the increased relativistic energy density due to 8N£jf + SN^f- [Both effects will 
be included in our numerical work]. It turns out that the mirror baryon acoustic 
oscillation effect is not only larger in magnitude, but has the distinctive feature of 
suppressing the higher odd peaks more than the even ones. 

Although the effect of additional relativistic energy density has been well studied 
in the literature (see e.g. ref.[26, 27]), and can be explored using existing CMB codes, 
the mirror baryon acoustic oscillation effect on the mirror dark matter perturbations 

4 The mirror photons can undergo diffusion (Silk damping). This would washout small scale 
inhomogeneities in the mirror radiation field just before mirror photon decoupling. However it 
should have very little effect on the ordinary CMB or matter power spectrum since p 7 > is a very 
small component to the overall energy density (assuming x ~ 0.5). Of course mirror photon diffusion 
would be expected to significantly damp the tail of the mirror CMB anisotropies, just like ordinary 
photon diffusion does for the ordinary CMB. 
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requires modifications. The relevant equations, though, are a straightforward gener- 
alization to the equations governing the perturbations of the ordinary baryons and 
photons. Our strategy is to numerically solve these equations, essentially using tech- 
niques developed in refs.[28] (see also [29, 30]). A very clear and helpful review is 
the one given by Dodelson[31]. 

Recall, the anisotropy spectrum today is characterized in terms of Cg. These 
quantities are the variance of the coefficients, ae m in the expansion of the photon 
temperature field in terms of spherical harmonics, Y^ m , i.e. 



The terms Cg can be related to the £ multipole moment, 0^, in the Legendre 
expansion of the Fourier transformed photon temperature field via the equation: 



where is the initial power spectrum of the metric perturbation with initial value 
$\ Finally the moments Oe(k, i] ) today can be related to the perturbations ©o(&, 
Vb(k,rj), $(/c,r/), ^/(k,rj), U(k,i]) near photon decoupling [where Vb(k,i]), ^/(k,rj), 
§(k,ri), Yl(k,r]) are the baryonic velocity, metric perturbations, and polarization 
tensor respectively]. The critical equation is[31]: 



where r is the optical depth for Thomson scattering and g(i]) = — fe~ T is the visibility 
function, which peaks near photon decoupling. The evolution of the quantities 0o, i>&, 

$ and II are governed by a set of linear equations, arising from the Boltzmann- 
Einstein equations. We assume standard adiabatic scalar initial conditions. The 
relevant equations are given in the appendix. 

It is important to note that the mirror dark matter model introduces only one 
additional parameter, x = T^/T 7 which is related to the fundamental Lagrangian 
kinetic mixing parameter, e via Eq.(2). The cosmological evolution of mirror dark 
matter, in the limit where x — > (i.e. e — > 0) exactly mimics cold dark matter. 
This is because mirror particles feel negligible pressure after the mirror photon de- 
coupling epoch, t' dec , and t' dec — > as x — > 0. As x increases from zero, differences 
begin to appear. Our job now is to determine what the observable differences are. 
To study these effects for the CMB one cannot simply choose a particular point 
for the parameters fl m h 2 , f2b/i 2 , h, ... from a fit assuming standard cold dark matter 




(5) 




(6) 
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and vary x. Doing this, for example, would modify the epoch of matter radiation 
equality, zeq + 1 = Q m /Q r , due to the additional contributions [Eq.(3)] to Q r . The 
matter radiation equality has been precisely constrained by the data and thus any 
modification to z E q by new physics needs to be compensated for by adjustments to 
the parameters (in this case, Q m h 2 ). In fact, what needs to be done is to examine 
parameter space where not only zeq is fixed, but also flbh 2 and 6 S (the angular size 
of the sound horizon at decoupling), since these quantities have also been precisely 
determined by the data. A similar situation has been noted when considering the 
effect of additional relativistic neutrino degrees of freedom[27, 26]. In this parameter 
space direction, the observable effects from varying x occur at small angular scales. 

It is reasonably straightforward to write a code to numerically solve the relevant 
set of equations to obtain the CMB anisotropy spectrum. For a given set of param- 
eters, flbh 2 ,Q m h 2 , h, comparison of our code with existing high accuracy codes, 
e.g. CMBFAST[32], confirms that our computation of the Ci values are accurate 
to within a few percent. This is sufficient for making a comparison of mirror dark 
matter with standard cold dark matter. 
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Figure 1: The anisotropy spectrum for mirror dark matter versus standard cold dark 
matter. The solid line is standard cold dark matter model with parameters described in 
the text (equivalent to mirror dark matter with x = 0), while mirror dark matter with 
x = 0.3 (dashed line), x = 0.5 (dotted line) and x = 0.7 (dashed-dotted line) are also 
shown. 
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In figure 1,2,3 we give our results for the CMB spectrum. We consider a flat 
Universe with the reference parameters fl m h 2 = 0.14, Q^h 2 = 0.022, = 1 — fl m , 
h = 0.70 [Q m = rib + fib']. These reference parameters are defined at x — 0. As 
discussed above, these parameters are adjusted as x is varied such that zeq, VL h h 2 and 
9 S are held fixed. [We also adjust the overall normalization by fixing the height of 
the first peak]. A scale invariant initial perturbation spectrum (Harrison-Zel'dovich 
and Peebles spectrum) is assumed and we have neglected reionization effect. Since 
we are interested in comparing mirror dark matter versus standard non-interacting 
cold dark matter (cosmologically equivalent to mirror dark matter with x — 0) small 
effects due to primordial tilt or reionization are not important to leading order. 
Figure 1 illustrates the expected agreement at large angular scales, as we vary x. In 
figure 2, we consider the small angular scale region of interest. In figure 3 we plot 
Fi(x) = C([x)jCi{x = 0) for several several values of x. 
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Figure 3: F?(x) = Cg(x)/Cg(x = 0) for x = 0.3 (dashed line) and x = 0.5 (dotted line), 
and x = 0.7 (dash-dotted line) are shown. 

Figures 2,3 clearly show the expected suppression of anisotropies at small angular 
scales, starting around the third peak. Interestingly, we see that the suppression is 
larger for the higher odd peaks than the even ones. These features can be readily 
understood. Odd peaks arise from compressions of the baryon-photon fluid, even 
peaks are rarefactions. When the gravitational driving force is suppressed, one ex- 
pects the odds peaks (the compressions) to be more affected than the even peaks 
(related effects occur when Q^h 2 is reduced). Furthermore, the differences only be- 
come apparent for the higher peaks because the suppression of power only occurs at 
small scales. 

Currently the most accurate measurement of the CMB damping tail has been 
made with the South Pole telescope [3]. These measurements show a slight damping, 
around ~ 2.5% at £ ~ 2000 c.f. predictions of the standard cold dark matter model. 
This damping provides an interesting hint that x ~ 0.4 [i.e. e « 2x 10 -9 from Eq.(2)]. 
In any case, these observations limit x ~ 0.5 [or e ~ 3 x 10~ 9 ]. It is anticipated that 
the PLANCK mission should improve the precision, which will probe e in the range: 
10- 9 £ e ~ 3 x 10~ 9 . 

In addition to CMB anisotropies the matter power spectrum can also be used 
to constrain parameters. However since small scales k ~ 0.1 h Mpc^ 1 have gone 
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nonlinear today, we consider the matter power spectrum on larger scales than this 
[linear regime]. It is straightforward to compute the power spectrum of matter, 

P(k) = 2n 2 5 2 *T*(k) (8) 
-"o 

where H = lOO/i km sec~ 1 Mpc~ 1 is the Hubble rate today and T(k) is the transfer 
function (see e.g. ref. [31] for details). In figure 4 we compare the obtained matter 
power spectrum for the various x values considered, for the same parameters used in 
figures 1-3. [Recall, Q m ,Qb and h are varied as x changes such that zeq, VL^h 2 and 
9 S are fixed]. As expected, deviations only occur on small scales as x increases from 
zero. This figure indicates that a rough bound of x ~ 0.3 — 0.5 could be extracted 
from galaxy surveys. Also note that similar results to our figure 4 have been obtained 
in the earlier study by Ciarcelluti[8]. 
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Figure 4: Power spectrum of matter for the same parameters as figure 1. As in figure 1, 
x = (solid line), x = 0.3 (dashed line), x = 0.5 (dotted line) and x = 0.7 (dashed-dotted 
line) . 



In conclusion, we have examined the implications of kinetically mixed mirror 
dark matter for CMB anisotropics. This dark matter candidate can potentially leave 
distinctive signatures on the CMB spectrum. We have found that the most important 
effects of kinetic mixing on CMB anisotropies is the suppression of the height of the 
third and higher odd peaks. This effect will be sensitively probed by the PLANCK 
mission in the near future. 
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Appendix - Linear perturbation theory with mirror dark matter 



The relevant equations governing the linear evolution of scalar perturbations in 
the Universe has a rich history starting with the work of Lifshitz in 1946 [33] and 
developed by many others, e.g. ref.[34]. For an up to date review see ref.[31]. As 
summarized in that review, the relevant equations governing the moments of the 
photon distribution (including the polarization field), which we consider numerically 
up to order t = 5 together with corresponding moments for neutrinos and baryonic 
matter perturbations, in the conformal Newtonian gauge, are: 
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(9) 



where f = —X e (l — Ypjubcrxa, Y p ~ 0.24 is the primordial helium mass fraction, 
a T is the Thomson cross-section and R = |^-. For the mirror sector, we have an 
analogous set of equations with 6^ — > Q' e , Qpe — > Q' pe N — > N' t (£ > 0), 5b — >■ 



f& ^ ^& and f,R—> f', i?'. Here r' = — X e /(1 — YL')n{/0"pa and i?' = |^-. Compared 
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with the standard cold dark matter model, the only additional parameter introduced 
is x = T^/T^ which is related to e via Eq.(2). [Recall these equations reduce to the 
equations governing standard cold dark matter when x — > and py — > p c ]. Finally, 
we have the two relevant Einstein equations: 



A; 2 ($ + 
k 2 $ + 3- ( $ - 



a 



-327rG'a 2 (p 7 2 + p<,iV 2 + p^ + Py© 2 ) 

47rGa 2 [p 6 <5 b + p b /5 b / + 4p 7 6 + 4p l/ iV + 4py6d + 4p l/ /^] 



(10) 
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where p Y = x 4 p 7 and p v = A e// (7/8)(4/ll) 4 / 3 p 7 , N eff = 3.046 + SN^ ff . For our 
application we can neglect N' e , p v > because we have negligible excitation of the mirror 
neutrino degrees of freedom. All derivatives in Eqs.(9,10) are with respect to confor- 
mal time, 77. The quantity, X e is the free electron fraction [X e = n e /riH where n# is 
the total number of hydrogen nuclei]. It obeys the Boltzmann equation[35, 31] 



ldX e 



a drj 



\l-X e )P-X 2 e (l-Y p )n b a( 2) 



C 



(11) 
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A a + A 27 + /3e 3 ^»/4T 7 • 



(12) 



Here eo = 13.6 eV is the binding energy of Hydrogen, A 27 = 8.227 sec and A a = 
//(3eo) 3 /[(87r) 2 (l — X e )nb(l — Y p )]. A similar set of equations will govern X e > (with 
Y p — > Y^, T 7 — > T' , lib — > rib'). Evidently, the latter depends on the primordial mirror 
helium mass fraction, Yp. This quantity can be computed solving the relevant mirror 
BBN equations, and for e ~ 10~ 9 , is [36] Y' p ps 0.85. [Note that Y' v is a slowly varying 
function of e if e ~ 10~ 9 ]. 

The equations must be supplemented with initial conditions. We consider the 
standard adiabatic scalar perturbations. We further assume that the initial pertur- 
bation of $ is drawn from a scale invariant Gaussian distribution with mean zero 
and variance, parameterized in the usual way: PJ, = (507r 2 /9/c 3 )<5^(f2 m /-Di(a = l)) 2 . 

The above set of equations, together with the Friedmann equation are numerically 
solved for k values on a logarithmically spaced grid between [fc min , k max ]. For our 
numerical work k m i n = 20/r^o and k max = 6000/^o- The Ci values are then obtained 
from Eqs.(6,7). 
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